Magnetic properties of an SU(4) spin-orbital chain 
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In this paper, we study the magnetic properties of the one-dimensional SU(4) spin-orbital model 
by solving its Bethe ansatz solution numerically. It is found that the magnetic properties of the 
system for the case of g t = 1.0 differs from that for the case of gt ~ 0.0. The magnetization curve 
and susceptibility are obtained for a system of 200 sites. For < gt < g s , the phase diagram 
depending on the magnetic field and the ratio of Lande factors, gt/g s , is obtained. Four phases with 
distinct magnetic properties are found. 



PACS number:75.30.Kz,75.10.Jm, 02.20.-a 

Recently, much attention has been focused on strongly 
correlated electrons with orbital degrees of freedom 
progress in experimental studies 
of transition metal and rare earth compounds such as 
LaMnC>3, CcBg and TmTe. Examples of spin-orbital 
systems in one dimension include quasi-one-dimensional 
tetrahis-dimethylamino-ethylene(TDAE)-C6o ||, artifi- 
cial quantum dot arrays [j)J and degenerate chains in 
Na 2 Ti 2 Sb 2 and Na 2 V 2 5 compounds 0|l]jl|. In 
these systems, the low-lying electron states have orbital 
degeneracy in addition to the spin degeneracy. This may 
result in various interesting properties associated with 
the orbital degrees of freedom in Mott insulating phases. 
For example, the magnetic ordering is influenced by the 
orbital structure which may change with pressure, or the 
magnetization is a nonlinear function of magnetic field 
even in the case of an isotropic exchange interaction of 
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The spin model with orbital degeneracy in a magnetic 
field was studied by means of effective field theories || 
recently. It has been shown that there exists critical be- 
havior under magnetic field. For a nonperturbative un- 
derstanding of the magnetic properties of the model, the 
Bethe-ansatz solution of the SU(4) model is applicable 
when the exterior magnetic field is applied since the to- 
tal z component of spin and orbital is a good quantum 
number. Although Ref. J(| studied the magnetic prop- 
erties in terms of the Bethe-ansatz method, their result 
involves only a special case because the Lande factor was 
not taken into account. In present paper, we study the 
model in the presence of an exterior field by solving the 
Bethe-ansatz equation numerically meanwhile taking ac- 
count of the Lande g factor. We obtain a rich phase 
diagram in comparison to what was obtained in || . Our 
results shed new light on the understanding of more real- 
istic systems. The quantum phase transition(QPT) 
concluded from such a system is speculated to be found 
in experiments. 

The one-dimensional quantum spin 1/2 system with 
twofold orbital degeneracy is modelled by 15 19] 
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where the coupling constant is set to unity. The Hamil- 
tonian (|l|) has not only SU(2)xSU(2) symmetry, but also 
an enlarged SU(4) symmetry Jig] . It was solved by the 
Bethe-ansatz approach, the obtained secular equation 
reads [fiT 



M 

2nI a = N0_ 1/2 (\ a )+Y,Qi(K-Xb) 

6=1 

M' 

+ E 6 '- 1 /2( A « - Mc), 
c=l 

M M' 

2irJ a = 22 6-1/2 (pa - Ab) + ^ 01 (/■*<» ~ Mc) 

6=1 
M" 

0-1/2 Oo - Vd) 

d—1 
I' 

2nK a = ^20_ 1/2 {v a - fi b ) + 6i (u a ~ v c ). 



c=l 



M" 



(2) 



6=1 



where 9 p (x) = —2 tan l (x/p). The quantum numbers 
{la, Ja, K a } specify a state in which there are N — M 
number of sites in the state |t), M - M' in |f), M' - M" 
in ||), and M" in ||). The A, /i and v are rapidities related 
to the three generators of the Cartan subalgebra of the 
SU(4) Lie algebra. The coupled transcendental equations 
(||) determine the energy spectrum of present model. To 
solve the Bethe ansatz equation numerically, one needs 
to take into account of the following propositions: 

• Given a N, Max(Af) = 3A^/4, Max(M') = N/2, 
Max(M") = N/4, and N - M > M - M' > 
M' - M" > M"; 

• {I a , J a , K a } are consecutive integers (or half inte- 
gers) arranged symmetrically around the zero, that 
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are integers or half integers depending on whether 
N-M-M',M-M'- M", M' - M" are odd or 
even respectively. 

The latter comes from the logarithm and the former is 
due to the property of the Young tableau. Once Eqs. (|^) 
are solved, the energy in the absence of exterior field is 
evaluated by 



M 

E (M,M' > M") = -J2 Tjr+x2 - 

a=l 1 X a 
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As the z component of total spin and that of to- 
tal orbital are given by S z = N/2 - M' and T z = 
N/2-M+M' -M" , the magnetization M z = g s S z +g t T z 
reads 

N 

M z = — {g s +9t) - Mg t - M (g s -g t )- M"g t , (4) 

where g s and gt are Lande g factors for spin and orbital, 
respectively. In the presence of the magnetic field H the 
energy becomes 



E{H, M, M', M") = E - HM Z . 



(5) 



In the absence of a magnetic field, the lowest-energy state 
is a SU(4) singlet where M" = N/4, M' = N/2 and 



M = 3N/2. In the presence of a magnetic field, how- 
ever, level crossing occurs and the state with M" — N/4, 
M' = N/2 and M = 3N/2 is no longer the ground state. 
Therefore the values M,M',M" for the lowest state de- 
pend on the magnitude of the applied magnetic field H. 
By analyzing the level crossing between the lowest-energy 
state and the others in terms of (5), we are able to cal- 
culate the magnetization by Eq. (4). 

We set g s = 2.0 and observe various situations for dif- 
ferent values of gt- The dependence of magnetization 
versus the magnetic field H is plotted in Fig. |l] for the 
cases of g t = 0.0 and g t = 1.0 respectively. Clearly, the 
magnetic properties of these two cases are quite different. 
We find that there exist three critical points in the mag- 
netization process, of which we will give interpretations 
on the base of symmetry analysis, the evolution of the 
Young tableau as well as the physics picture by means of 
the mean-field analysis later on. 

For an easier understanding, we first give an explana- 
tion by means of the Young tableau for 8-site system. 
In Table I, we present the whole 15 different multiplets, 
their energy in the absence of a magnetic field, and the 
corresponding magnitudes of magnetization. Each mul- 
tiplet corresponds to a distinct Young tableau. 



TABLE I. The ground state energy Eq in vanished exterior field, the z-component of the total spin S z , that of the total 
orbital T z , and magnetization M z of all possible configuration of an eight-site system 



Multiplet 


[2 4 ] 


[3221] 


[3311] 


[332] 


[4211] 


[422] 


[431] 


[44] 


[51 3 ] 


[521] 


[53] 


[511] 


[62] 


[71] 


[8] 


Eo 


-14.92 


-13.73 


-13.84 


-13.41 


-12.76 


-12.85 


-12.04 


-11.30 


-10.83 


-9.95 


-10.25 


-7.41 


-7.99 


-4.00 





T z 





1 





1 


1 


2 


1 





2 


2 


1 


3 


2 


3 


4 


s z 





1 


2 


2 


2 


2 


3 


4 


2 


3 


4 


3 


4 


4 


4 


M z (g t = 0) 





2 


4 


4 


4 


4 


6 


8 


4 


6 


8 


6 


8 


8 


8 


M z (g t = 1) 





3 


4 


5 


5 


6 


7 


8 


6 


8 


9 


9 


10 


11 


12 



The application of an exterior field will make each en- 
ergy level corresponding to a multiplet that carries out a 
SU(4) irreducible representation split into Zeemann sub- 
levels. If g t = 0.0, the energy level of an SU(4) singlet 
labelled by the Young tableau [2 4 ] will cross with the 
sublevel of the highest weight state labelled by the Young 
tableau [3311] when the applied exterior field reaches the 
value of 0.27. If the exterior field increases further this 
state becomes the ground state. The magnitude of mag- 
netization plays the role of the slope of energy curves 
determined by Eq.(|). When H = 0.635 its energy level 
will cross with that of the state labelled by Young tableau 
[44]. For H > 0.635 the spins are totally polarized and 
the magnetization is saturated. 







ff>0.27 































ff>0.635 



When g t — g s = 1, the level crossing occurs to the high- 
est weight states labelled by Young tableaux [3311] , [332] , 
[422], [44], [53], [62], [71], and [8] for the 8-site system. 

It is clear that the exterior field only breaks the spin 
SU(2) symmetry and the orbital SU(2) symmetry re- 
mains when g t = 0. The magnetization process is similar 
to that of traditional Hciscnbcrg model but the inter- 
mediate state are those labelled by the Young tableau 
[mi, m2, m3, TO4] with mi = m,2 and 7713 = m,4, i.e., 
M' = 2M", N + M' = 2M. This implies that the evo- 
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lution of Young tableau caused by the exterior field is 
realized by moving a couple of boxes lying in the third 
and the fourth row to the first and the second row. Unlike 
in the transitional Heisenberg model where an increasing 
exterior field flips the down spins to up spins one by one, 
it always flips a pair of parallel spins that constitutes an 
orbital singlet. As a result, present spin-orbital model 
presents larger magnetic susceptibility meanwhile it un- 
dergoes two phases (i.e., phase I and III in Fig. ||) versus 
the exterior field. 

Because the orbital SU(2) symmetry is also broken 
for nonvanishing g t , the magnetization process becomes 
complicated when 0.0 < gt < 2.0. The whole process 
undergoes four regimes as long as gt < g s - The first 
regime consists of strong spin polarization and weak or- 
bital polarization (SSP-WOP) processes, where the level 
crossing occurs to the states labelled by four-row Young 
tableaux [mi, m%, m^, m.4]. The contributions to M z of 
both spin and orbital are positive in this regime. The sec- 
ond regime corresponds to the level crossing between the 
state labelled by three-row Young tableaux, i.e., m 4 = 0. 
Because the contribution of spin to M z is positive, but 
that of orbital to M z is negative, we call it the regime 
of strong spin polarization and orbital anti-polarization 
(SSP-OAP) regime. The third regime is related to the 
level crossing occurring among the states labelled by two- 
row Young tableaux (7713 = 7714 = 0) where spin degree of 
freedom is frozen and orbital degree of freedom takes part 
in the polarization process only. We call it the regime 
of spin frozen and orbital polarization (SF-OP). Eventu- 
ally, the magnetization reaches saturation and then both 
the spin and orbital become frozen (SF-OF) which is the 
phase IV indicated in Fig. |2|. For g t = 1.0 the three criti- 
cal values of the magnetic field are — 0.31, iJf = 0.93, 
and HI = 4.0. 

Define R = gt/g s , we plot the phase diagram of H 
and R in Fig. g. The critical lines that separate the 
aforementioned regimes are plotted in Fig.(^). Both the 
spin and orbital degree of freedom are frozen in the phase 
IV. There is a simple relation for the boundary line be- 
tween SF-OP and SF-OF regimes, R = 2/H, which arises 
from the competition between the state related to the 
Young tableau [N - 1,1] and [N]. The state in SSP- 
WP regime are quantum liquid of spin-orbital mixture 
while SF-OP is orbital quantum liquid. If R = 0.0, 
the state in SSP-WOP has SU(2)<g>U(l) symmetry, while 
that in SF-OP has SU(2) symmetry @. When R ap- 
proaches to unit, however, the third phase SF-OP dis- 
appears, the phase SSP-OAP directly transits to phase 
SF-OF at point H c = 2.0 then it undergoes only three 
distinct regimes which will be explained again in the fol- 
lowing. 

In the above we interpreted the phases with the help 
of the properties of Young tableaux. It can also be un- 
derstood on the basis of mean-field analysis. In the pres- 
ence of magnetic field, at first, the Zeemann term in the 



Hamiltonian brings about three processes: orbital flip- 
ping up, spin flipping up or both spin and orbital flip- 
ping up simultaneously, which is the phase I indicated in 
Fig. 0. When the magnetic field exceeds certain value 
for < g t < g s , the three flipping processes, ||) — > ||), 

I J.) — * |t) and ||) — » |t) are stopped. In phase II, the fol- 
lowing flipping processes take place, |J.) — > |t), |t) ~ > It) 
and ||) — > |t). The last one is a process of spin flipping 
together with orbital anti-flipping, which suppress the 
susceptibility (slop of the magnetization curve) in phase 

II in comparison to that in phase I. As is known |l^] that 
the ground state of the present model in the absence of 
exterior field is an SU(4) singlet whose physics picture 
can be considered as the simultaneous singlet of the spin 
and orbital, as well as the product of spin and orbital. 
Starting from the Neel state of the product of spin and 
orbital, the mean-field analysis tell us that a pair of par- 
allel spins will favor their orbital to be anti-parallel, or 
vice versa. This makes it easy for us to understand the 
spin flipping but orbital anti-flipping process in phase II. 
Increasing the magnetic field will make it cross to the 
edge of phase II where the spins are completely polar- 
ized. Thus in phase III, the orbital start to flip back 

t) — » |t)- Once all orbital become completely polarized 
it goes into phase IV. 

If g t = 0, the phase I transits to phase III directly. 
Because the orbital SU(2) is not broken, orbital up and 
down are equivalent and always half were in orbital up 
and half in orbital down among the spin flipping pro- 
cess in phase I. Increasing the magnetic field until all 
spins are completely polarized, the orbital liquid state 
still stay in a liquid state and it will not flip even if in- 
creasing the magnetic field further. So it will not go into 
phase IV. If gt — g s , the pure spin flipping and orbital 
flipping have the same probability in phase I. In phase II 
the spin flipping with orbital anti-flipping dose not occur 
and pure spin flipping and orbital flipping have still the 
same probability. Thus spin polarization itself will not 
happen unless spin and orbital are completely polarized 
simultaneously. Therefore it transits from phase II into 
phase IV directly. 

In summary, we have investigated the magnetic prop- 
erties of SU(4) spin-orbital chains by numerical solutions 
of their Bethe-ansatz equations. If the exterior field is 
decoupled to the orbital degree of freedom, i.e. gt = 0, 
we found that there exists a critical point of the quan- 
tum phase transition which separates it into two regimes. 
If the magnetic contribution of the orbital is not ne- 
glected, however, there will be three critical points and 
four regimes that demonstrate different magnetic prop- 
erties. Though our calculation is made in one dimension, 
the analysis based on the symmetry and symmetry break- 
ing has no restrictions in dimension. For a high dimen- 
sional systems (two or three dimensional) , it is speculated 
to have similar features. 
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FIG. 1. The magnetization of SU(4) spin-orbital chain versus external 
to a* — 1 .0 and the dashed line corresnonds to o* = 0.0. 
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FIG. 2. Phase diagram of magnetic field, H and the ratio of Lande g 
(SSP-WOP), II. (SSP-OAP), III. (SF-OP), and IV. (SF-OF) in the regio: 
figure corresponds to g t = 1.0. 
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